Intrinsic parton transverse momentum in next-to-leading-order pion production 
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We study pion production in proton-proton collisions within a pQCD-improved parton model in 
next-to- leading order augmented by intrinsic transverse momentum (fcjj of the partons. We find 
the introduction of intrinsic transverse momentum necessary to reproduce the experimental data in 
the CERN SPS to RHIC energy range, and we study its influence on the so-called K factor, the 
ratio of the NLO cross section to the Born term. A strong pr dependence is seen, especially in the 
3 — 6 GeV transverse momentum region of the outgoing pion, where nuclear effects (e.g. the Cronin 
effect) play an important role. 

PACS numbers: 24.85.+p,13.85.Ni,25.75.Dw 
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I. INTRODUCTION 

Today's collider facilities raise the interest in testing 
perturbative QCD (pQCD) at work, and in searching for 
phenonema beyond its capability. However - at least 
for the experimentally available transverse momentum 
region - pQCD parton models underestimate thepro- 
duction of mesons in proton-proton (pp) collisions [J, Q , 
even at next-to- leading order (NLO) In order to re- 
store the consistency with the data, two methods were 
proposed at leading order (LO): inclusion of the intrinsic 
transverse momentum of partons 00, and/or an effec- 
tive correction factor (K factor), accounting for higher 
order contributions 0, 0. The physical background is 
to account for both, the missing higher order perturba- 
tive corrections and radiation effects. The latter is not 
present in electron-proton processes, however, it becomes 
important in pp collisions [6j . 

In this paper we present the first results on pion pro- 
duction in pp collisions applying a NLO pQCD parton 
model with intrinsic transverse momentum, taken from 
a Gaussian distribution with width (kj_). In Section II 
the intrinsic transverse momentum is introduced into the 
formalism, and the appropriate NLO expressions are pre- 
sented. The necessity of such an extension is demon- 
strated at y/s = 27.4 GeV in Section III. Next, we dis- 
play the best fit values of the width of the intrinsic trans- 
verse momentum distribution at NLO level for available 
pp — > 7r + X experiments, in the energy range 20 GeV 
< y/s < 200 GeV. Similarly to Ref. % we study the 
pion production in the 2 GeV < pr < 6 GeV transverse 
momentum region, where nuclear effects are considered 
to be important. 
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Finally, we extract the ratio of the NLO cross section 
to the Born term (K factor) at different energies and 
transverse momenta, and study its dependence on the 
amount of the included intrinsic transverse momentum. 
This way we provide a numerical foundation to the cor- 
rection factors used in LO calculation 0. We demon- 
strate, that a leading order calculation with a fitted, 
cm. energy, transverse momentum and scale dependent 
K factor and additional intrinsic transverse momentum 
reproduces well the full NLO results at higher energies 
and momenta, and can be used as a fast method to get 
a reasonable estimate of a full NLO calculation. 



II. MODEL 

In order to extend the applicability of the original, infi- 
nite momentum frame parton model to smaller trans- 
verse momenta, we introduce the intrinsic transverse mo- 
mentum of the partons |fj. We write the four- momenta 
of the interacting partons (a and b) as @ 



Pa \%a ~r /— , K_L,ai %a 



K ±,a 



2 2x a y/s 

k 2 



2 2x a y/s 



Pb 



, s ' s K X,b r \fs , k l,b 
[x b — + - A;_l )6 , ~Xb— + - 7= 

2 ZXby/S 2 2Xby/S 



In this notation x, the momentum fraction carried by the 
parton, becomes a parameter. The apparent fraction is 
x — kj_/(xs), however, for practical applications at high 
energy (p T > 3 GeV; yfs > 40 GeV; {k\) < 2 GeV 2 ), the 
distinction has a negligible (< 5%) effect. Furthermore, 
we require that the longitudinal direction of the partons 
does not change sign due to the transverse momentum, 
i.e. x > k±/y/s. 

The starting point of our calculation is factorization: 
the hadronic cross sections up to a power correction may 
be written as a convolution over hard partonic (pQCD) 
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processes, 
da 



dyd p T 



/dz 
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dcr 



fa/pi x a,Q; k± t a)fb/p(%b, Q\ fc_L,fe) — J Ar/c (z c , Qj) , 

where dcr /dt is the partonic cross section of the reaction 
a + b — > c + d (LO with condition <5(1 + (t + u)/s)), or 
a + 6 — > c+<i+e (NLO with fixed z c ) and (at fixed scales) 
is the function of the partonic Mandelstam variables only. 
In order to avoid singularities due to the intrinsic trans- 
verse momentum we use regularization 
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u - M 2 /2 , (3) 



with M = 1.8 GeV. The factorization is done at the fac- 
torization scale Q, where the parton content of the initial 
proton is determined by the parton distribution function 
fa/p (PDF). For simplicity, we assume a Gaussian depen- 
dence of the PDFs on the intrinsic transverse momentum, 
with a width (k'jj, 
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We note, that such a separation may be viewed as a first 
approximation to the unintegrated PDF, where recent 
studies indeed found a shape close to Gaussian [ljj. Fi- 
nally, the hadrons are created collinearly with the outgo- 
ing parton c with momentum fraction z c at fragmentation 
scale Qf. The partonic cross section explicitly depends 
on three scales, Q, Qf and the renormalization scale Q r . 

In principle, the scales can be determined such that 
the final result has the minimum sensitivity to them ||, 
and usually are set to be equal. However, in this paper 
we use the results of a previous study | ll| of pp and 
pA data, where the reproduction the Cronin effect in pA 
reactions ^3 imposed such scales, that the corresponding 
(fci)«2GeV 2 . 

We note that there is some ambiguity in the choice 
of scales. In the literature typical scales are fixed to 
hadronic or partonic variables, kpt or kpt/z c , respec- 
tively, where k is an 0(1) number. Other choices are also 
possible, e.g. an invariant scale, Q 2 — k 2 siu/ (s 2 +P+u 2 ) 
as proposed in Ref. Q. However, in our case we found 
that this choice is equivalent to kpt/z c . 

At NLO level and no intrinsic transverse momen- 
tum Eq. J2J) is usually rewritten with variable change 
(x a ,Xb,z c ) — > (v,w,z c ), where t = —(1 — v)s and u — 
—vws, as 
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with the proper kinematical boundaries, J being the 
Jacobian of the transformation (1/J = u(l — v)w for 
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FIG. 1: Comparison of experimental data fl^l to the NLO 
pQCD parton model result in p + p — > ir + + X reaction at 
Eiab = 400 GeV, with and without intrinsic transverse mo- 
mentum. 



(fcl)=0), and da NLO /dv is the sum of 2 
cross sections 13], 
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J Born (f) \ 

— — — S(l-w) + ^±iic« b ^(s, i, w, Q, Q r ,Q f ) ,(6) 

dv 7T 

with renormalization scale Q r (chosen to be equal to the 
factorization scale Q). In this paper, however, we use 
Eq. J2J directly, since for (k 2 ^) ^ the momentum frac- 
tion x cannot be expressed analytically from v and w, 
while the inverse transformation can be done. 

Several codes are available for calculating jet cross sec- 
tions at NLO level [l3L lh"H . Here we have chosen to ex- 
tend the one by the Aversa group calculating the par- 
tonic cross sections at next-to-leading-log level, with the 
intrinsic transverse momentum distribution. The calcu- 
lations presented in the following sections were performed 
with the MRST-cg PDF [l| and KKP FF [ig parame- 
terization. 



III. RESULTS 

A. Comparison to data 

First, we demonstrate the importance of the intrinsic 
transverse momentum to reproduce experimental data in 
the transverse momentum region 3 GeV < y/s < 6 GeV. 
In this part, we use the scales proposed in Ref. to 
make a direct comparison to the 400 GeV FNAL exper- 
iment 17]. Fig. 1. shows, that at y/s = 27 A GeV the 
NLO calculations of pion production in pp collision un- 
derpredict the experimental data by a factor of 2 using 
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reactions from 19.5 to 200 GeV. Different symbols refer to 
different experiments, see [Dl l20jl. The solid line represents 
the average value of a LO calculation 



the scale parameters Q — Q r = Qf = pr/2 and ne- 
glecting intrinsic transverse momentum (as in 0). De- 
creasing the scales (but still keeping them hard for the 
applicability of pQCD) the agreement can be improved, 
however, without intrinsic transverse momentum, even if 
lower scales are chosen, the data are still underestimated. 

The lower line in Fig. 1. presents the NLO calculation 
with a partonic intrinsic transverse momentum distribu- 
tion of width (fc^) =1 GeV 2 , and shows a nice agreement 
with the data. We note, however, that there is a delicate 
interplay between the choice of the scales and the intrin- 
sic transverse momentum (k\) : needed to reproduce the 
data 01- Typically, increasing the scales increases the 
value of (fc 2 ). 



B. Intrinsic transverse momentum 

In this section we summarize the results on the in- 
trinsic transverse momentum width of the partons in 
the nucleon fitting the NLO pQCD calculations to the 
data, generalizing the LO scale choice of 

[12 to Q = 

Q r = kpt/z c , Qf = upt, where n is a 0(1) number. 
In this work, we fix k = 2/3 at NLO level. Our previ- 
ous study [TjJ also showed, that the reproduction of the 
Cronin peak in pA collision requires the width of the in- 
trinsic transverse momentum distribution to be on the 
order of (k\) w 2 GeV 2 at energies yfs ~ 30 GeV, which 
is achieved by the above choice of scales. A similar value 
of (kj_) was extracted from the experimental analysis of 
jet-angle distribution [l9| . 



Analyzing the pp — > n + X experimental data [rH Il9l 
|20| . we deduced the best fit value (kj_) for each exper- 
iment, similarly to what is shown in Fig. 1, minimizing 
the x 2 (-C /T — 1) (data over theory) ratio in the range 
3 — 6 GeV. The result is presented in Fig. separately 
indicating the runs from different experiments and shows 
a need for a considerable amount of partonic transverse 
momentum. We also checked, that similarly to the LO 
results pj, the extracted width does not depend on the 
charge of the pion. 

Usually, changing the order of a calculation requires 
changing of the underlying scales. Comparing NLO re- 
sult to the previous LO ones pj one notices that in order 
to keep the average transverse momentum width, we had 
to increase the scales,dictated by the experimentally ob- 
tained Cronin peak [ll|. Keeping LO scales (k = 1/2) 
would lead to a substantial reduction of the width, origi- 
nating in the mechanism of NLO graphs to automatically 
generate transverse momenta. It is remarkable, that fit- 
ting to the nuclear reaction data requires the same width 
independent of the order of the calculation used! 

Recent dAu data at RHIC [2l| also indicate that at 
\fs = 200 GeV more transverse momentum is necessary, 
than (fc 2 ) w - 0.5 GeV 2 , indicated in Fig. 2. This 
shows, that the scale parameter k possibly may also de- 
pend on the energy, and a higher, k — 4/3 scale with 
(fc^) = 2.5 GeV 2 reproduces well both the pp and the 
dAu data without jet quenching effects j2^. The depen- 
dence of the scale parameter on the cm. energy and its 
consequences will be studied in a separate paper. 
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FIG. 4: Ratio of the pionic K factor at (k\) = 1 GeV 2 (left) 
v/i = 20 GeV to 1800 GeV. 



C. K factor 

Since most calculations (especially, for nucleus-nucleus 
collisions) are still based on LO, it is useful to provide a 
well-founded K factor for these faster calculations. Fig.|3| 
shows the ratio of the full NLO calculation to the Born 
term with no intrinsic transverse momentum. Indeed, 
at high energy and transverse momenta K n approaches 
the well known value of 2, however, in the low transverse 
momentum region it has a strong px dependence. As 
a first approximation, for yfs > 60 GeV, the pionic K 
factor can be taken energy independent [23j |. 

Since the intrinsic transverse momentum may have a 
different effect on the Born term than on the higher or- 
der (HO) processes, the latter contributing the dominant 
part, it is worth to study the dependence of the K factor 
on the width (kj_). We present this behavior in Fig. 0] 
indicating the ratio R = K({k^))/K(0) at {k\) = 1 and 
2 GeV 2 (with k — 2/3 scale parameter) and at different 
energies from -y/s = 20 to 1800 GeV. While at low ener- 
gies the ratio R shows a large decrease with increasing 
(kj_) (mainly due to the efficiency of the intrinsic trans- 
verse momentum at low energies in enhancing the Born 
term faster than the HO corrections), from ^/s > 60 GeV 
its value is always above 1, indicating, that HO correc- 
tions are gaining more from transverse momentum, than 
the Born term. The maximal enhancement is peaked at 
j)i«2-3 GeV, with 10% correction at (kj) = 1 GeV 2 
and 20% correction at = 2 GeV 2 , while for p T > 

4 GeV, the transverse momentum driven correction van- 
ishes, validating LO calculations in this range with (kj) 
independent K factors |^.ITl|. 
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and 2 GeV 2 (right) to K factor at {k\) = at energies from 



IV. CONCLUSIONS 

In this paper we introduced initial transverse momen- 
tum distributions into the parton-based description of 
hadron production in NLO level pQCD calculations and 
demonstrated, that such an extension is necessary even 
in NLO to reproduce the pp experimental data. From 
the analysis of most experiments in the 3 GeV < pr < 6 
window a width of intrinsic transverse momentum distri- 
bution on the order of {k\) ss 2 GeV 2 was fitted. The 
obtained precision of the description of pion production 
in pp collisions is high enough to find possible collective 
effects in nuclear collisions. 

We presented the pionic JsT-factor (full NLO cross sec- 
tion to the Born term) for several energies and transverse 
momentum values, and found a pronounced dependence 
on pt in the 3 GeV < pr < 6 window, where nuclear 
effects show up most prominently. Furthermore, we in- 
vestigated in detail the modifications of the Born and 
higher order terms due to the presence of an intrinsic 
transverse momentum and concluded, that above certain 
energy in the above mentioned pt window the higher or- 
der contribution raises faster, than the Born term and 
hence the i^-factor increases compared to its value with- 
out intrinsic transverse momentum. For pt > 4 GeV 
this enhancement may be neglected, justifying (k^_) inde- 
pendent _ftT-factor calculations, however, at smaller trans- 
verse momenta the correction goes up to 20%, just in the 
order of the measured nuclear enhancement. The de- 
tailed study of the nuclear enhancement, and its role in 
fixing the scales of an NLO calculation will be carried 
out separately. 
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